ASYMPTOTIC BEHAVIOR OF BRANCHING 
POPULATIONS BEFORE EXTINCTION 
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Abstract. Under the assumption that the initial population size 
of a Galton- Watson branching process increases to infinity, the 
paper studies asymptotic behavior of the population size before 
extinction. More specifically, we establish asymptotic properties of 
the conditional moments (which are exactly defined in the paper) . 



1. Introduction and the main result 

We consider a Galton- Watson branching process {X n } n > , 

(1-1) Xn+l = y^^n.j, 

3=1 

where X n denotes the number of offspring in the nth generation for a 
population starting from K offsprings, X = K, and throughout the 
paper the initial size of population K is assumed to be a large value. 
Such a type of branching process can be a model of real population of 
animals, insects etc., and the main results of our study can have ap- 
plications to analysis of real populations arising in biology (e.g. Jagers 
[TT] . Haccou, Jagers and Vatutin [10], Jagers and Klebaner [13]). For 
other study of branching processes with a large initial population size 
see also Borovkov pQ, Klebaner [IT], Klebaner and Liptser [T8] . 
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The study of branching populations before extinction has been ini- 
tiated by Jagers p2] and then resulted in papers of Jagers, Klebaner 
and Sagitov [3] and [15]. The approach of these papers [H] and [15] 
is based on analytic techniques for studying paths to extinction with 
the following analysis of asymptotic behavior of these paths. 

The present paper presents an alternative way to study asymptotic 
behavior of large populations before extinction, and the approach of 
the present paper is based on diffusion approximations of the original 
branching process with large initial population as well as a series of aux- 
iliary processes. Those diffusion approximations are then used to study 
asymptotic behavior of conditional moments of a population size before 
extinction as it is explained below. The approach of our paper (includ- 
ing diffusion approximations, asymptotic expansions and sample path 
techniques) remains correct for much wider classes of branching pro- 
cesses than that traditional branching process and includes for instance 
bisexual Galton- Watson branching processes [I] , [5] and different type 
of controlled ^-branching Galton- Watson processes (which need not 
be subcritical, as it is assumed in the paper later). The ^-branching 
processes have been introduced by Sevastyanov and Zubkov [21J, and 
intensively studied in many papers (e.g. Bruss [2], [3], Gonzalez, Molina 
and Del Puerto [7j, [8], [9], Yanev [22], Zubkov [23]). 

The main results of the present paper are presented in Theorem 11.11 
below. 

Assume that £ n j, (n=l,2,. . . ; j—1,2,. . . ) have the same distribution 
for all n and j and are mutually independent, and there exists the 
second moment E££ • < oo. Denoting m = E£ n j and (5 2 = var(£ nj ), 
assume that m < 1. Under this last assumption the extinction time of 
the branching process always exists with probability 1. Let r = tk be 
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that moment of extinction. The random variable tk is a stopping time 
associated with the sequence {X n } n > . We assume that the family of 
all stopping times {tk} (for different values K) is defined on a filtered 
probability space {fl, & Q , F = (& 0iK ), P}, ^o,a- C i^x+i C . . . C #q. 
(The meaning of the index will be clear later.) 

The paper studies asymptotic behavior of the branching population 
before extinction as K increases to infinity, and the main result of our 
study, formulated below, as well as the analysis of the paper, use the 
notation x for asymptotic equivalence between two main parts of ex- 
pansion. The notation is used in order to reduce irrelevant background 
explanations and to avoid multiple using of lim in different senses or 
expansions with remainder, where it is not significant. For example, 
relations (11.2)1 and ( II. 3p (see formulation of the theorem below) should 
be read as follows: For any sufficiently small positive e and 5 there 
exists a large integer K such that 

P {(1 - 5)X[ U2TKl Em l{luiTKi - lu2TKi) 

< E(X[ UlTKj | X [u2TKl ) < (1 + 8)X[ U2TKj Em l ^i-^^} > 1 - e, 
P {(1 - 5)X[ UiTK] Em l{lu2TKl - luiTKi) 

< E(X[ U2Tk] I X [U1TK] ) < (1 + tf)X[ wJ Em^^J-L^J)} > 1 - e , 

and 

P {(1 - 5)K l m lluiTKl 

< E(X[ U1TK} I tk) < (1 + 5)K l m 1 ^ } > 1 - e. 

In the places where it is required and looks more profitable (e.g. 
Section EJ), the explicit form of asymptotic expansion with remainder 
is used nevertheless. 
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Theorem 1.1. Let < u\ < u 2 < 1 be two real numbers. Then, as 
K — > oo, 



1.2) 



E{X[ UiTk1 | V^x^^Em^^J-W), 
E{X[ U2TkJ I ^jlx^jEm^Hw^ 

and 

(1.3) E{X[ UiTifJ | r.JxifW^, 



where \_z\ is the notation for the integer part of z. As K — > oo, j v 
converges in probability to the constant c = — j^^. 

The proof of the main result is based on the following lemma. 

Lemma 1.2. For any finite- dimensional vector {X^, X i2 , . . . , X in } , 
1 < Zi < i 2 < ■ ■ ■ < i n < oo, 

„(X h - Km h X i2 - Km' 2 X % - Km in 

i im p 1 — r- n? - Xl> r- n? - X2 ' • • • ' — r~ 7W - Xn 



= P{0 H <xi,0i 2 <x 2 ,...,8 in < x n }, 

where {0i,9 2 , ■ ■ ■} is a Gaussian sequence with E9j = and cov(9j, 
6 j+n ) = m n var(6' i )+nm J+n - 1 ; var(# i+1 ) = m 2 var(^)+m J ', var(^) = 1. 

Lemma [1.21 is known from the literature, and its proof can be found 
in Klebaner and Nerman [TH] . For the purpose of the present paper we, 
however, need in an alternative proof of this lemma, which follows from 
the asymptotic expansions presented here. Furthermore, the proof of 
Theorem 11.11 requires the intermediate asymptotic expansions obtained 
in the proof of Lemma 11.21 rather than the statement of Lemma 1.2 
itself. 
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In this paper, simple asymptotic representations for all conditional 
moments before extinction are obtained. The most significant conse- 
quence of this analysis is a so-called invariance property of the condi- 
tional expectations. This property is discussed in Section [61 

The main idea of the method is as follows. The random sequence 
{X n } n > is approximated by appropriate random sequences 

n>0 

(a < 1), as a tends to zero. For each fixed a we define stopping times 
T a ^ K (for different values K) associated with the process Yn ^ . T~a,K is 
assumed to be measurable with respect to the cx-field ^ a ,K C where 
&a — ^K>i^a,K, and & a C J^o- F° r that fixed a the sequence r aK 
converges (in definite sense) to £(a) as K — > oo (the details are given 
in the paper). Then knowledge of the behavior of YutI, k , < u < 1, 
for which we have the corresponding relationship, enables us to study 
the behavior of its limit as a tends to zero. This limit is just X UTK , 
< u < 1. Other assumptions associated with definition of X n and 
that of the associated processes X { n\ Y^ a) and other processes are given 
in the next section. 

The rest of the paper is organized as follows. In Section [2] we intro- 
duce the auxiliary stochastic sequences X^ and Yn^ and the stopping 
times associated with these sequences. The elementary properties of 
these random objects are studied. In Section [3] we continue to study 
the properties of the sequences Xffl and Yn a \ Specifically, it is shown 
that these sequences are upper and lower bounds for the branching 
process X n , and these bounds are tight as a — > 0. These properties 
are then used in order to prove the convergence results in the next 
sections. In Section H] we derive asymptotic expansions and prove the 
convergence lemma to the Gaussian process, the parameters of which 
are explicitly defined in the formulation of Lemma [L2l In Section [5] we 
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prove Theorem 11.11 Last Section [6] discusses application of the main 
results of this study and establishes the invariance property. 



2. Stopping times and auxiliary processes associated with 
the Galton- Watson process 

In this section we approach the stopping time tk, the extinction mo- 
ment, by introducing a parametric family of stopping times {r aj ^}, 
depending on the two parameters a and K. Specifically, for any real a, 
< a < 1 and integer K 

(2.1) r a , K = mi{l : X t < [aK\}, 

where \_aK\ is the integer part of aK. The stopping time T a ,x as well 
as the associated with these parameters a and K other corresponding 
random variables defined below are assumed to be measurable with 
respect & a ,K C Tk, and for two different values a\ and a 2 , < a 2 < 
ai < 1, we have JP ait K C ^ a2 ,K- If a < 1 is fixed and K ly K 2 are 
distinct, K\ < K 2 < 00, then we have & a ,K\ C ^ a ,K 2 - Then the 
two-parametric family of u-fields {J^ a ,K} is increasing in the following 
sense. For any < a 2 < «i < 1 and integer K\ < K 2 < 00 we have 

a\,K\ — ^ a 2 ,K 2 - 

In accordance with this family of stopping times (12.11) . consider a 
family of processes X^ a K = Xj satisfying the recurrence relation (for 
notational convenience the additional index K is not provided): 

{ X W \ 
[aK\,J2j n A, X^=K. 

The processes Xj K are assumed to be adapted with respect to the 
ex-fields & a ,K- m addition, the processes X^ a K are assumed to be 
measurable with respect to the the wider a-field J^o- Specifically, if 
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there are two processes X^ and with different a\ and say 

< a 2 < Oi < 1, then both of these processes X^g and -ST.- j£- are mea- 
surable with respect to the cr-field ^o,K, and, of course, with respect 
to the cr-field ^ a2 ,K- All of these processes with different a are defined 
due to representation (12. 2p . This means that the processes X^ K are 
actually defined after their stopping times as well. For different a\ and 
02 (0 < a 2 < a\ < 1) the processes Xjj) and X^ are 'coupled' until 
the stopping time r ai ,K, i-e. until that time instant their sample paths 
coincide, but after the time instant r ai K these processes are decou- 
pled i.e. their paths become different. But the coupling arguments 
can be used nevertheless: after the time instant r ai ^ with the aid of 
Kalmykov's theorem [16] we have X^ > st X^ , j > r ait K (see the 
next section for details). 

Some mathematical details about these processes can be found in the 
next section. The similar coupling arguments hold for the processes 
YjK defined later, which are derivative from the processes X^ (the 
further details can be found in the next section). 

Let us transform ( 12. 21) by adding and subtracting the term [a^FTj . 
To this end we use the following elementary property of numbers: 
max{a, b} — a = max{0,6 — a}. Also there is used the fact that 
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Xn > [aK\ for any n. Then, we have 
fx^-laK^ + laK] 



x 



(a) 



max | [aK\ , £ nJ J> - [aK\ ] + [aK\ 
(2-3) f 4^ L^J 

= max ^ 0, Znj " E fnj + (l - CnjjJ > + L^J 



max 



.(a) _ v (a) 



Hence, denoting Yn = Xn — from ( 12. 3ft we obtain 



max 



F n ( J = max <J 0, £ E (' ' 

j=[aK\+l j=l 

[o*J 



i=i i=i 



\J'=1 3=1 

\J'=1 3=1 

where £^ • = ^j+^j £^ 2 >--- are independent and identically 

distributed random variables having the same distribution as £ n> j), and 
lt } = I { n\K) = I{ Y,f=l Cn,j > EK J (1 - L tj )} is the notation used 
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9 



Thus we have the new family of processes Y„ , which is assumed, as 
mentioned before, to be measurable with respect to J^o an d given on 
the same probability space {Q, J^o, P}- Recall that a stopping time 
T a ^K and the sequence 



are assumed to be adapted with respect to the cr-field ^ a ,Ki and the 
family of these cr-fields {^ a ,K} is increasing when a decreases and K 
increases. 

It is known that as K — > oo, ^ converges to m n in probability (see 
Klebaner and Nerman [19]). Using this result it is not difficult to prove 
that, as K — > oo, 



(2.5) 





converges to max{0,m n — a] in probability, 



i 



K 



(a) 



n 



converges to max{0,m n — a} 



K 



I)?\K) converges to Xn+i = Xn+i(m, a) 




in probability, 



as well as, 



(iv) 



T at K converges in probability to £(a) = min{/ : m < a}. 



The proof of (i) is postponed to the end of Section [31 The proofs of 
(ii) — (iv) are similar to the proof of (i). 
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3. Properties of the sequences Xn and Y^ 

The study of this section we start from the properties of the random 
vectors (12.51) . Let ai, a,2 be two numbers, and < 02 < a\ < 1 and K 
is fixed. Then, in the suitable probability space for all events u> G Q 
and n > 

(3.1) Y^\u)<Y^\u). 
Indeed, consider two random vectors 

(3.2) {y^\y^\..) 
and 

(3.3) {y^\y^\..) 

Consider the stopping times {r ai ,K, ^ ai ,x} and {r a2t K, ^o, 2 ,k} associ- 
ated with the sequences (13. 2j) and (13. 3j) . Since for fixed K, JP aij K f= 
^ a2 ,K, then T auK {u)) < T a2iK (u). 

According to the definition of the sequence X„ (see (12.21) ). on the 
suitable probability space containing ^ a2 ,K we have the correspondence 

(3.4) X^\u>)=X^\u), • = 1,2,..., 

and at this stopping time r ai ^ K we have X^ k (uj) > Xt^ k (uj), and 
therefore according to Kalmykov's theorem [Id] : 

A i d.st A i j * — TolK) T oi,iC + - 1 ! • • • j T a<z,K ) ■ ■ ■ j T 0,K- 

Therefore, in a suitable probability space 

(3.5) Xl ai \uj)>X\ a2 \u), i = T auK ,T auK + l,...,T a2 ,K,---,T Q ,K- 

Thus, we showed 

X- ai \u) > X\ a2 \u), i = 1,2,.. .,r 0l ,K,r 0li K + l,...,T a2 , K , ■ ■■,r Q>K - 
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From this correspondence (13.41) and ( 13.51) according to the definition 
of the sequence Yn (see (12.41) ) on the same probability space we have 

(3.6) Y^\u) + [ ai K\ = Y^iu) + [a 2 K\, i = 1, . . .,r ai , K - 1, 

and therefore up to time r ai ,K — 1 the inequality Y^ ax \uj) < Y^ a2 \uj) is 
obvious. At time instant r ai ^ we have Yt^ k {uj) = 0, while Y^ K (u>) is 
nonnegative in general. The further behavior of the processes Y^\u>) 
and Y^ a2 \u) after time T aii # is specified by coupling arguments, where 
the initial inequality Y^ ai \u) < Y^ a2 \u) before the stopping time r aij K 
remains true after this stopping time as well. If for some i = io, 
Y^Ioj) = Y^ 2 \uj) (= 0), then the both processes are coupled un- 
til ii > i . If after time ii, Y^?\(u) becomes positive, then we again 
arrive at the inequality Y^\(u) < Y^_\(u), and so on. 

Taking into account that according to the definition Xn coincides 
with X n , we obtain the inequality 

(3.7) y(«)M<i n M<ii fl) M, 

being correct for all u G £1 and all n > as well as for any initial 
population K and any a. This inequality is also tight as a — > 0, because 
according to the definition of the above sequences, Yn ^(w) — X n ^(to>) 
for all n. 

Let us now prove the above properties (i) — (iv). Find the limit in 
probability of -gr- as K — > oo. Notice first, that according to (J272J) — ^- 
converges to max{a, m} in probability, and according to Wald's equa- 
tion [6], p. 384, ^ converge to the same limit max{a, m}. Therefore, 

y( a ) 

converges to max{a,m} — a = max{0,m — a} in probability, and 
^ converges to max{0,m — a}. Now, assuming that for some k it 
is already proved that -j^- converges to max{0,m — a} in probability 
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and £ converges to max{0, m — a}, by induction we have as follows. 
If m < a then -k^- converges to in probability and -h converges 
to 0, and consequently, 

[ 3=1 3=1 

The last is true because 



^ «; f 



mEY^ 



K K ^^' J K 

3=1 3=1 

]7;y( a ) y( a ) 

Therefore, according to (j2.4p — vanishes, and vanishes in 
probability. Therefore, the assumption m k < a is not the case. Hence, 
assuming that converges to m — a in probability, where m > a, 
we have the following: 

EY {a) 

lim — -§±i = max{0, (m fc - a)m - a(l - m)} 

fe+1 



max{0,m + - a}. 



■i^(a) 

Thus, as X — > oo, converges to max{0, m n — a} in probability, and 
(i) is proved. Notice, that (ii), (Hi) and (w) follow together with (i). 
All these claims are closely related, and their proof is similar. 

Notice also, that the convergence of -gr- to max{0, m™ — a] in prob- 
ability means that in a suitable probability space, the sequence Yn K ^ 
converges almost surely to max{0,m n — a}. 



4. Asymptotic expansions and the proof of Lemma 11.21 



Pathwise inequalities (13.71) and < hold for any initial size K and 
any a. Therefore the appropriate normalized sequences -g^- and — 
converge to the same limit in probability as K — > oo. If there exists 
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the limit in distribution of Xn ~^ in as K — > oo, then because of the 
equality y ~ o) ^E y ~ o) = xP-ex™ ; anc i tne inequality F„ (a) (a;) < X n (u) 
for all a > and all u G f2 (see ref. ( 13.71) ). there are also the limits 

,. ., . „ y(»)_I?v('') _ \-( a )_Pv( a ) -r r , „ 

m distribution of " ^= " and " ^ n as K — > oo and a — > 
independently. That is, one can let — > oo before a — ■> 0, or converse. 
Notice, that the limiting distribution of Xn ^ Xn has been obtained in 
[19] . and it also follows from asymptotic expansions obtained in this 
section. 

It follows from the results of Section [3] that, as K — > oo, t^k con- 
verges in probability to 



£(a) = min{/ : m l < a}. 



and hence, in the case where K increases to infinity first, £(a) = 
P-limx-»oo T~a,K (P"lim denotes a limit in probability). It is known 
(see e.g. Pakes [20]), that converges to the constant c = — in 
probability. This result of Pakes [20] can be proved by different ways. 
The advantage of the proof given below is that it remains true for 
more general models than the usual Galton- Watson branching process, 
resulting in the justice of the results of the paper for general models 
as well. For instance, one can reckon that a bisexual Galton- Watson 
branching process starting with K mating units is considered, where 
m now has the meaning of the average reproduction mean per mating 
unit (see Bruss [1]). For the relevant result related to the ^-branching 
processes see Bruss [3], Theorem 1. 
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For large Xn = K we have as follows: 

t k : = inf{* G N : X t = 0} 
= inf{t e N : X t < 1} 

=inf|tGN:^<l 
X A 



inf it 6 N : Yl 



X n 1 

n=l 



Now note that, as K —> oo, each fraction * n converge to m in prob- 
ability. Indeed, 

X n X„ K 



(4.2) 



X n -i K X n _i 



According to |19j . — > m ?l in probability as K — > oo. Therefore, the 
fraction (14. 2[) converges to = m in probability for any n. 

On the other hand, by virtue of Wald's identity [6], p. 384 we obtain: 

^ EX n _ E 1 Cn-ij _ mEX n _! _ ^ 

EX n _i EX n _! EX„„ 1 

So, according to (14.21) and (14. 3p . the limit in probability of the fraction 

/" as K — > oo and the fraction of the corresponding expectations 



are the same. 



From (14.31) we therefore obtain: 

(4 . 4 ) B (£)-^*-n 



EX TO 

X n I K *-\ EX n _! 

n=l 



m*. 



lim KEm TK = 1. 



So, 



Similarly to (14.41) . we also have that yj^J converges to m* in probability 
as K — > oo for any integer t. So, from (14.11) and (14. 4p we have the 
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similar limit as in the case of expectations, i.e. for any positive e and S 
there exists integer K large enough such that P{|ZTm TA ' — 1| > 5} < e, 



i.e. Km TK — > 1 in probability as K — > oo, and thus converges to 
— logm * n probability. 

From (12.41) for Yj^, j = 0, 1, . . ., we obtain the following equations: 
(4.5) 



;=i 6\/& ] K 



i=l 



where • oo is assumed to be 0. 

Assuming that K increases to infinity, and dividing both sides of 
(14.51) by large parameter &\fK we have the following expansions 



(4.6) 



y j+1 - mXj+lIj 



6, 



evK 

y( a ) 



m 



Xj+i 



K 



\aK\ 



i=l 



or 



(4.7) 



Xj+i ■ 



Xj+i 



Y {a) - mY {a) 



V' 



6 



\ 



1 -m 



For j — 0,1, ... , 1(d) — 2, £(a) > 2, one can remove the term Xj+i from 
the both sides of (147]) . 
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Therefore, for j = 0,1,..., £(a) — 2, the left-hand side of (14. 7p can 
be transformed as follows: 

Y (a) - mY {a) 
r- +1 mi 



Y}1\ - mEY} a) + mEYj a) - mlf } 



&VK 

Y}Z\ - EY}°\ yW - Elf > aVK / 
— m— ——[ 1 - m 



where {9j } is a Gaussian sequence. (The values of the parameters of 
this Gaussian sequence are not discussed here.) 

In turn, for j = 0, 1, . . . ,£(a) — 2, £(a) > 2, the right-hand side of 



(14.71) is transformed as 



1=1 v 2 = 1 



(4.9) 



e,/y/ a) 



0V^--^(l-m 

where {C,} is a sequence of independent standard normally distributed 
random variables. 

Therefore for j = 0, 1, . . . , £(a) — 2 from (14. 8 p and (14.91) we have: 



(4.io) Ci-t<'x0Vy 



where yj°^ = 



(a) 
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The analysis of (14.101) is standard. According to the definition ~E9j = 
0. Therefore, rewriting (14.101) as 



(4.H) ^-mflj a) +0VK! 8) . 
we obtain: 

(4.12) E(^\) 2 = jimE(m^+C/ 



' Y ( a ) 



K 



Taking into account that E ^CjyVj^J = ^ , y l f^ ~> m ^ - a, as X -> oo, 
we obtain 

(4.13) var(^.+ ) 1 ) = m 2 var(#j a) ) + m 3 - a. 

Next, from (14. 10|) we have: 

cov(0f ,0^) =E6f6% 1 



lim E I mefef + QOfJ^ 

K-+oo \ 33 3 3 V K 



mvar(^j a ' ) ) + m-' — a, 



and it is easy to show by induction 

cov(0< a \ O = m" var(^) + £ m^Ey't- 



i=l 

&y j+n -i — tn - a, 

where j ' + n < £(a) — 1. Assuming now that a — ► 0, we obtain the 
convergence of the sequence 

Xi - Km X 2 - Km 2 
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to the Gaussian process {9±, 9 2 , . . .} with mean and covariance func- 
tion 

cov(0 i; 6 j+n ) = m n vai(dj) + nm J,+n-1 , 

(4.14) 

var(^ + i) = m 2 var(^) + m j , var(#i) = I. 
This implies the statement of Lemma II. 21 

5. Proof of Theorem 11.11 

Let us now study equation (14.101) more carefully. Let U\ and u 2 be two 
real numbers, < U\ < u 2 < 1- Assume that K is so large that the 
probability P{\r at K — £(a, K)\ > e} is negligible (e > is an arbitrary 
fixed value, K is large enough), where £(a, K) is a (not random) integer 
number. Such a number does always exist for any given a since, as 
K — > oo, r a t x converges to £(a) in probability. 

For large K we have the following two expansions: 



where j/f^ T , in the right-hand side of equations (15.11) and (15. 2p . 



i = 1,2, are correspondingly replaced by % l? ,K . It is worth noting as 
follows. Relations (15.11) and (15.21) are written in the form of an asymp- 
totic expansion. The left-hand sides of these expansions are Gaussian 
martingale-differences, while the right-hand sides are the expressions 
with large parameter K. Since the probability P{|r a A - — £(a, K)\ > e} 
is negligible (e > is an arbitrary fixed value, K is large enough), the 
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expansion with the given right-hand side is correct. From ( 15. lft and 
(15.21) we obtain as follows: 

and for any continuous function /(•) 



(5.3) 



For example, from ( 15.31) we obtain: 



(5.4) 



Now estimate the conditional expectation E <j [Y^ T k \ 
For brevity let us introduce a random vector 



Y, 



(a) 



\.U2T a ,K\ 



We have 



E1Y 



r (a) 



eIe(y, 



"2,r ai K 



(5.5) 



,(a) 

V u i T a,K J 



E 



(EfK (a) 



T a ,K 



T~a,K 



Ul,U2,T at K 



EY^ ., i = l,2. 



The last equality of the right-hand side of ( 15.51) is a consequence of 
conditional independence of Y^ t , and Z Ul)U2iTojJf , that is for any 
given event {r a> K — k}, the random variable Y$ k , and random vector 
Z uliU2i fc are independent. (I5.5P holds true also in the case of a = that 
will be discussed later. 
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Next, using the notation tk = tq^k let us prove that 

(5.6) cov {0\u lTK \+i - ttt0|u irK j, 0[uaTK-J+l ~~ m ®[u 2 T K \) — >■ 
as K — ► oo. 

Notice, first (see relation (14.141) ) that cov(9j,9j +n ) vanishes as n — > 
oo. Consequently, by the total expectation formula, 

(5.7) cov(#|„ 1TK j , 0[_ uirjf j +n ) = E (cov(0|_ UlTir j , 9\u 1TK j +n \ tk)) 

vanishes as n — > oo, where here in relation (15. 7p and later the notation 
for cov(9i UlTK j , 9i UlTK \ +n | Tk) or another similar notation means the 
conditional covariance. Taking into account that, as K — ► oo, tk 
increases to infinity in probability and u 2 — U\ > 0, the difference 
|_ u 2Tk"J — L^i'^J increases to infinity in probability too. Hence, by 
virtue of (15. 7\i one can conclude that cov(#|_ UlTK j , 9\_ U2TK j) vanishes as 
K — > oo. Therefore, as K — > oo, ~Ei9\, aiTK \9\u 2TK \ is asymptotically 
equal to E0|_ UiriC jE0|_'u 2riC j, anc l (15-61) follows. In addition to (15. 5ft and 
(15.61) we have also the following. Since the sequence {Q} consists of 
independent standard normally distributed random variables, then as 
K — ► oo 

(5-8) cov (C[ ul r K i, C[u 2 r K \) -> 0. 

This is because cov (Cluit^j , C[u 2TK i I r K ) = \{[u 1 t k \ = [u 2 r K \}, and 
the last vanishes in probability as K — > oo. 

Assuming that a vanishes we need a stronger assumption than above. 
Specifically, we assume that K is so large that the probability 

T a ,K ~ £{a, K) 



lo S X 



> e 



log m 

is negligible for all < a < clq (e > is an arbitrary fixed value, K is 
large enough), where a$ < 1 is some fixed small number. Such a large 
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number K does always exist, since as K — > oo and a vanishing, 
converges to — j^-^ in probability. Then, letting a — > in f)5.4p in view 



(a) 



of pathwise inequalities (13.71) and Zj— < ^ we have 



^Ur K J (^L«2T K J+l - ™9[u 2 t k \) Cf Ul 



(5.9) 



2 ' '21 

U 2 T K \ ■ 



Taking into account (15.51) . (15. 6p and (15.81) and conditional indepen- 
dency of X [uiTKi , (9[ ulTK j + i - m6 [uiTKi ) and (,\_ U2TK \, and passing to 
the appropriate conditional expectations, from (15. 9p we obtain: 



E\X 



i 

\_UlTK J 



X 



[U 2 T K \ 



(5.10) 



E 



X 



ElX 



i 

\U\TK J 



X [u 2 t k \ , TK 



21 



E {0[ Ul r K \+i - tng|_ UlTK j) 

E (9[u 2 T K \+l - ^L"2T K j) 



Thus, to this end our task is to determine the asymptotic of 



for large Returning to basic equations (15. ip and (15. 2p . we have 



2/ 



/ 



E 



C[UiT aiK \ 



Y, 



(a) 



2; 



L«iT"a,ifJ 



V 



(5.11) 



E < E 



(a) 



2/ 



~K 



T~a,K 



E 



where (i UlTaK \, i = 1,2, are standard normally distributed random 
variables. As a vanishes, from (15.111) we obtain 



(5.12) E (9 [UiTKl+1 - m6 lUtTK] ) 21 x Em<^E(C Lu 



TI(\ , 



\2l 
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Therefore, (15.101) can be rewritten 

This proves the first equation of (11. 2p . The proof of the second equation 
of (11.21) is similar. 

Consider basic equation (15.11) again, rewriting it as follows: 

Assuming that as a vanishes we have: 

(5.14) X\,a 1TK \Clu 1 r K \ x K (^L«itkJ+i ~ m ^L«i^/d) 2 • 

Therefore, taking into account that X\_ UlTK \ and C[u 1 t k \ are condition- 
ally independent, from (15.141) we obtain: 



E ^ X [uiT K \ 



( 11.31) is proved. 





( IU1T K \+1 ~ m6 luiTKi ) 21 




E (([uiT a , K \ 


r K ) 21 



6. Discussion 

The aim of this section is to present the main results in convenient 
form for application to analysis of real populations. In this section we 
also establish a so-called invariance property. 

Let, when K is large, e be a relatively small (positive or negative) 
parameter having the following meaning. The population size at time 
[u 2 tk\ is assumed to be equal to [(1 + e)Km U2tK \ , t K = — 

The meaning of this value is the following. The factor Km^ U2tK ^ 
is the expected size of the population at time L^TkJ > an d the factor 
1 + e represents a parameter of relative deviation from the expected 
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population at that time moment. Then, from Theorem 11.11 we obtain, 
that for large K 



(6.1) 



logE(X[ ttirxJ | X [u2TKl = [(l + e )Km^\) 
= l\og(K + Ke) + lu x t K logm + o(l) 
= l\og(K + Ke) -l Ul hgK + o(l). 



In real computations the term \og(K + Ke) can be replaced by e + \ogK 
if e is sufficiently small. 

The result similar to (16.11) can be obtained for the conditional ex- 
pectation of (ll.3p . Specifically, for large K write 

\ogK 



(6.2) 



1 + e 



logm 



con- 



(16.21) has the following meaning. As K — > oo, the fraction - o ^ K 
verges in probability to — j^^, and therefore, as K is large, the factor 
1+e is a parameter for relative deviation from the expected value of 
extinction time. Then, 



IokE^X? 



(6.3) 



- [U1T K ] 



tk = ~ 



logiT 
logm 



/ \og(K + Ke) - lux log K + o(l) 



As we can see the right-hand sides of (16.11) and H 6 . 3 1) coincide. That is 
for any given relative deviation 1+e the asymptotic conditional expec- 
tations are invariant. 
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